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I. 


SUMMARY 


A.  INTRODUCTION 

While  a  great  deal  of  work  has  been  done  on  laser  effects 

with  pulsed  infrared  lasers,  ^  only  recently  have  UV  lasers  of 

(2) 

suitable  energy  levels  been  developed.  In  the  IR,  the  inter¬ 

action  between  a  pulsed  laser  and  an  aluminum  surface  in  air  is 
dominated  by  the  formation  of  a  plasma,  which  totally  absorbs 
the  laser  energy.  The  plasma  then  re-radiates  a  large  fraction 
of  this  energy  in  the  UV.  Since  aluminum  is  a  very  good  reflector 
in  the  IR,  but  a  good  absorber  in  the  UV,  this  leads  to  higher 
surface  coupling  than  would  be  the  case  if  the  plasma  were  absent. 
In  addition,  the  plasma  generates  a  significant  impulse  level  on 
the  target  by  remaining  at  high  pressure  for  a  significant  period 
of  time.  With  a  UV  laser,  where  the  photons  are  much  more  ener¬ 
getic  and  where  the  basic  absorptivity  of  the  surface  is  higher, 
the  phenomenology  may  be  different. 

Under  IRAD  funding,  AERL  measured  the  thermal  coupling  of 
KrF  (249  nm)  and  XeF  (353  nm)  laser  radiation  to  metal  targets, 

1.  McKay,  J.A. ,  Schriempf,  J.T.,  Cronburg,  T.L.,  Eninger,  J.E., 
and  Woodroffe,  J.A.,  "Pulsed  CC>2  Laser  Interaction  with  a 
Metal  Surface  at  Oblique  Incidence,"  Appl .  Phys.  Lett.  3_6 
125  (1980). 

2.  Jacob,  J.H.,  Hsia,  J.C.,  Mangano,  J.A.  and  Ronki,  M.J.,  Appl. 
Phys.  50,  5130  (1979) . 
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(3  4 ) 

as  well  as  the  impulse  coupling.  '  These  results  indicated 
that  damage  by  UV  lasers  should  be  very  efficient.  Similar 
results  have  been  obtained  by  other  groups. 

All  of  the  experiments  to  date  have  been  with  relatively 
low  pulse  energy  lasers.  As  a  result,  measurements  at  high 
fluence  have  been  carried  out  with  small  spots.  In  extrapolat¬ 
ing  to  a  large  spot  situation,  it  is  important  to  determine  what 
are  the  dominant  physical  effects.  The  spot-size  scaling  will 
be  very  different,  for  instance,  if  the  phenomena  at  high  fluence 
are  purely  vapor  dominated,  as  opposed  to  plasma  dominated. 

This  report  covers  work  performed  under  a  DARPA-funded 
prooram  for  Analysis  of  UV/Visible  Laser  Effects  Data.  The  pur¬ 
pose  of  the  program  was  to  determine  what  effects  are  important 

2 

at  high  fluence  levels  lOJ/cm  ),  to  help  determine  spot  size 
scaling.  This  knowledge  is  required  for  extrapolation  to  system 
level  energies  and  for  proper  design  of  damage  experiments  on 
larger  machines.  Specifically,  the  program  focused  on  two  key 
aspects:  a)  transient  vaporization  of  aluminum  into  vacuum; 

and  b)  plasma  formation  in  aluminum  vapor  at  the  temperature  and 
density  predicted  by  (a) . 

3.  Woodroffe,  J.A.,  Hsia,  J.  and  Ballantyne,  A.,  "Thermal  and 
Impulse  Coupling  to  an  Aluminum  Surface  by  a  Pulsed  KrF  Laser," 
Appl.  Phys.  Lett.  3j6  ,  14  (1980). 

4.  Duzy,  C. ,  Woodroffe,  J.A.,  Hsia,  J.  and  Ballantyne,  A., 
"Interaction  of  a  Pulsed  XeF  Laser  with  an  Aluminum  Surface," 
Appl.  Phys.  Lett.  3_7 ,  542  (1980). 

5.  Rosen,  D.I.,  Mittledorf,  J.,  Kothandaraman ,  G. ,  Pirri,  A.N., 
and  Pugh,  E.R.,  AIAA  Paper  No.  AIAA-80-1321 ,  presented  at  13th 
Fluid  Plasma  Dynamics  Conference,  Snowmass,  Colorado,  July  1980. 

6.  Walters,  C.T.,  "Ultraviolet  Laser  Effects  Experiment,"  Report 
No.  BCL-G7029,  Battelle  Columbus  Laboratories,  Columbus,  Ohio 

(1980)  . 
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This  report  is  divided  into  three  sections  and  an  appendix. 
Section  I  presents  the  rationale  for  the  choice  of  tasks  (a)  and 
(b) )  and  contains  a  summary  description  of  the  conclusions  of  the 


work  described  in  the  other  two  sections.  Section  II  describes 
a  calculation  for  vaporization  of  aluminum  into  a  vacuum  under 
conditions  where  steady  state  may  or  may  not  be  reached.  Section 
III  presents  a  calculation  for  plasma  formation  in  aluminum  vapor. 
Finally,  Appendix  A  contains  a  description  of  the  AERL  IRAD  work 
on  KrF  and  XeF  laser  effects. 

B.  RATIONALE 

2 

At  low  incident  fluence  levels  (<1  J/cm  ) ,  the  interaction 
between  the  UV  laser  light  and  the  surface  results  simply  in 
warming  up  of  the  surface  (see  Appendix  A) .  As  the  fluence  in¬ 
creases,  the  phenomenology  changes.  First,  isolated  surface 
defects  may  reach  the  vaporization  point,  resulting  in  the  appear¬ 
ance  of  impulse.  Next  the  surface  melts  in  bulk  during  the  pulse. 
At  still  higher  levels  (~10-30  J/cm  )  the  surface  reaches  the 
vaporization  temperature  during  the  pulse.  At  this  point  sig¬ 
nificant  impulse  levels  appear.  The  impulse  is  expressed  in 

terms  of  the  coupling  coefficient  C,  which  represents  the  amount 

2 

of  impulse  intensity  (dyn-sec/cm  )  obtained  per  unit  of  fluence 
(J/cm  ) .  At  the  same  time,  the  thermal  coupling  (expressed  in 
terms  of  a,  the  ratio  of  energy  left  in  the  target  after  the 
pulse  to  the  incident  energy)  begins  to  go  down  because  the 
vapor  carries  away  with  it  a  great  deal  of  energy.  The  vaporiza¬ 
tion  process  is  transient,  since  the  surface  is  heating  up  through¬ 
out  the  pulse. 
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2 

Finally,  at  very  high  fluence  (>>30  J/cm  )  two  phenomena 
are  dominant:  strong  (near  steady-state)  vaporization,  and 
plasma  formation.  A  plasma  may  form  before  or  after  strong 
vaporization  begins.  In  the  strong  vaporization  regime,  high 
values  of  C  (~3-10  dyn-sec/J)  are  observed  (see  Appendix  A) .  If 
no  plasma  is  formed,  then  the  effect  of  spot  size  scaling  is 
small.  Because  the  temperatures  are  low  and  therefore  the  sound 
speeds  are  low,  information  about  the  existence  of  spot  edges 
propagates  slowly  to  the  center  of  the  spot  (speed  ~ 1.5x10  cm/ 
sec) .  Rather  small  spots  (<1  cm  area  for  a  pulse  length  '  2y 
sec)  are  suitable  for  near  one-dimensional  experiments.  The  re¬ 
sults  can  then  be  easily  extrapolated  to  the  fully  one-dimensional 
situation  at  operational  energies. 

If  a  vapor  plasma  is  formed,  however,  then  it  can  result 
in  values  of  C  much  greater  than  10  dyn-sec/J,  and  can  also  re¬ 
sult  in  higher  thermal  coupling  to  the  surface  (as  has  been 
observed  in  the  IR,  for  both  impulse^  and  thermal  coupling  ^ 8  ^  )  . 
In  the  IR,  the  plasma  medium  is  ambient  air.  In  the  UV,  however, 
the  principal  interest  is  in  vacuum  interaction.  Plasma  forma¬ 
tion  and  behavior  are  then  much  more  complex,  since  the  medium 
(vapor)  has  varying  initial  density  and  temperature.  If  a 
plasma  is  formed,  the  temperatures  and  sound  speeds  are  much 

7.  Reilly,  J.,  Ballantyne,  A.,  Woodroffe,  J.A.,  "Modelling  of 
Momentum  Transfer  to  a  Surface  by  Laser  Supported  Absorption 
Waves,"  AIAA  J.  \1_,  1098  (1979). 

8.  Pirri,  A.N.,  Root,  R.C.and  Wu,  P.K.S.,  "Plasma  Energy  Trans¬ 
fer  to  Metal  Surfaces  Irradiated  by  Pulsed  Lasers,"  AIAA  J. 

16,  1296  (1978). 
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greater,  resulting  in  an  order  of  magnitude  bigger  spot  size 
requirements  for  near  one-dimensional  experiments. 


Given  the  interest  in  extrapolating  to  very  high  fluences 
and  large  spots,  we  chose  to  perform  analyses  on  (a)  the  tran¬ 
sient  vaporization  process  to  determine  impulse  levels  if  only 
vaporization  is  important,  and  to  obtain  initial  vapor  condi¬ 
tions  as  a  function  of  incident  flux;  and  (b)  plasma  formation 
time  in  aluminum  vapor,  for  the  initial  conditions  obtained 
from  the  vaporization  calculations.  With  the  latter  informa¬ 
tion,  we  can  then  determine  whether  a  plasma  will  be  present  or 
not  at  a  particular  flux.  A  task  that  remains  for  future  work 
is  the  integration  of  the  two  analyses  to  obtain  a  predictive 
capability  for  the  surface/vapor/plasma  system  for  both  thermal 
and  impulse  coupling. 

C.  SUMMARY  OF  CONCLUSIONS 

The  two  analyses  are  described  fully  in  Sections  II  and 
III.  Only  the  conclusions  will  be  given  here. 

1 .  Section  II,  Transient  Vaporization 

1)  Generally,  heat  conduction  into  the  bulk  metal  and 
vapor  efflux  represent  equally  important  mechanisms 

for  removing  absorbed  power  from  the  surface  or  vapor 

interf ace . 
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2)  The  surface  heat  up  time  varies  inversely  with  the 
square  of  the  absorbed  power  and  becomes  small  com¬ 
pared  to  the  incident  pulse  duration  at  incident  XeF 
fluxes  >100  MW/cm^. 

3)  The  impulse  coupling  coefficient  C  is  one  parameter 

relevant  to  understanding  damage  mechanisms  associated 

with  pulsed  devices  interacting  with  a  surface.  At 

100  MW/cm2  absorbed  and  1  nsec,  C  ~3  for  a  -20%  (XeF) 

and  ~6  for  a  -40%  (KrF) ,  corresponding  to  an  incident 

2  2 

fluence  of  500  J/cm  and  250  J/cm  ,  respectively. 

2 .  Section  III,  Plasma  Formation 

1)  Plasma  formation  was  not  predicted  at  an  incident  flux 

2 

of  25  MW/cm  for  any  of  the  sets  of  rate  constants 

2 

considered  in  this  work,  while  at  100  MW/cm  plasma 
was  formed  very  rapidly  regardless  of  the  assumptions 
used  in  determining  the  rate  constants.  Between  the 
two  extremes,  the  relative  magnitude  of  the  photo¬ 
ionization  and  electron  impact  ionization  rates  de¬ 
termined  whether  a  plasma  is  present. 

2)  Experiments  with  1-2  usee  pulses  at  high  fluences 
(>100  J/cm  )  will  be  in  the  plasma  regime  and  there¬ 
fore  large  spots  and  a  laser  with  large  pulse  energy 
are  necessary  to  do  the  required  one-dimensional 
experiments . 
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3)  Several  improvements  can  be  made  in  the  model,  such 


as:  inclusion  of  a  non-Maxwellian  electron  distribu¬ 

tion  for  determining  the  electron  impact  rate  constants; 
a  better  estimate  for  electron  impact  ionization;  a 
more  detailed  calculation  for  the  photo-ionization 
cross  section  from  the  4p  state  of  aluminum. 

4)  Re-radiation  from  the  plasma  while  probably  a  small 
effect  needs  to  be  included  in  the  model  in  future 
work. 
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II.  TRANSIENT  VAPORIZATION  FROM  A  SURFACE  INTO  VACUUM 


A.  INTRODUCTION 

It  is  possible  to  direct  energy  into  a  surface  at  high 
power  density  (e.g.,  >  1  MW/cm  )  and  thereby  induce  rapid  surface 
vaporization.  High  power  lasers  represent  one  such  energy  source. 
There  have  been  several  studies  of  the  ensuing  processes  when 
there  is  a  surrounding  air  environment  '  '  '  ' .  The  case  of  a 
surface  initially  surrounded  by  a  hard  vacuum  is  less  well 
studied  and  is  the  object  of  this  section.  The  external  energy 
source  is  taken  to  be  pulsed  and  absorbed  only  at  the  surface 
(or  phase  interface  with  the  vapor  flow  region) .  Absorption  in 
the  vapor  phase  is  not  considered  even  though  it  can  be  important 
in  processes  such  as  laser-induced  breakdown.^ ' That  repre¬ 
sents  a  separate  process  which  will  be  pursued  in  Section  III. 

The  basic  methodology  to  be  outlined  is  fairly  general, 
but  it  is  convenient  to  choose  a  particular  material  in  order  to 
render  the  discussion  more  concrete.  An  aluminum  surface  will  be 
considered  herein.  The  primary  mechanisms  for  removing  the  ab¬ 
sorbed  power  from  the  surface  are  typically  heat  conduction  into 
the  dense  phase  behind  the  surface  and  efflux  of  vapor  from  the 

9.  Shui,  V.,  Rival,  B.  and  Weyl,  G. ,  "Effect  of  Vapor  Plasma  on 

the  Coupling  of  Laser  Radiation  with  Aluminum  Targets,"  J. 
Quantum  Spectroscopy  and  Radiative  Transfer  20,  627,  (1978). 

10.  Weyl,  G.,  Pirri,  A.,  Root,  R. ,  "Laser  Iqnition  of  Plasma 

Off  Aluminum  Surfaces,"  AIAA-80-131 9,  13th  Fluid  and  Plasma  Dy 
namics  Conference,  July  1980. 
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surface,  carrying  with  it  the  heat  of  vaporization.  Theoreti¬ 
cal  modeling  of  these  two  processes  will  be  discussed  in  turn. 
Then  the  discussion  will  proceed  to  coupling  them  with  the  trans¬ 
ient  vapor  flow  away  from  the  surface. 

B .  MODEL 

2 

At  the  absorbed  power  fluxes  of  interest  (>  1  MW/cm  ),  the 
aluminum  surface  will  very  quickly  heat  up  to  a  temperature  well 
above  the  melting  point.  Thus,  both  solid  and  liquid  phases  will 
generally  be  involved  in  the  real  process,  with  a  moving  phase 
interface  between  them.  This  has  been  avoided  in  the  analyses 
done  thus  far  by  imagining  that  there  is  only  a  single  dense 
phase,  a  liquid.  The  heat  of  fusion  for  aluminum  is  small  com¬ 
pared  to  the  heat  of  vaporization.  It  has  also  been  assumed 
that  the  thermal  properties  are  constant.  This  is  decidedly  not 
true  as  seen  from  Figure  1.  The  approximations  just  noted  rep¬ 
resent  considerable  simplications ;  however,  it  can  be  expected 
that  the  general  trends  will  be  representative.  Heat  transfer 
into  the  dense  phase  is  modeled  by  the  linear  heat  conduction 
equation: 

9T  _  .  5 2T 

PlCSL  3t  »  2  ^ 

3  x 

where  p ^  is  the  liquid  density,  c ^  is  its  specific  heat,  and 
is  its  thermal  conductivity.  The  possibility  of  convection 
in  the  liquid  (driven  by  vapor  recoil,  surface  tension,  etc.^11^) 
will  be  neglected. 

11.  Palmer,  H. ,  "The  Hydrodynamic  Stability  of  Rapidly  Evap¬ 
orating  Liquids  at  Reduced  Pressure,"  J.  Fluid  Mechanics 
75,  487  (1976) . 
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The  surface  will  recede  as  the  aluminum  boils  off.  This 

is  depicted  in  the  x-t  diagram  in  Figure  2,  where  the  dense  phase 

zone  lies  to  the  left  and  vapor  flows  in  the  positive  x-direc- 

tion.  It  is  convenient  to  employ  a  moving  reference  frame 

attached  to  the  surface  as  this  simplifies  imposing  the  liquid/ 

vapor  phase  jump  conditions.  To  this  end  introduce  n  =  -Jx  +  x_(t)] 

s 

with  x  =  -x  (t) at  the  surface.  The  heat  conduction  equation 


becomes 


3T  _  .  3T  ,  3  T 

3t  *s3n  ^  2 

3n 


where  =  k^/p^c^  is  the  thermal  diffusivity  and  ^s  =  dxg/dt. 
Initially  T  =T^  for  all  n  > 0  and  also  this  provides  the  appro¬ 
priate  boundary  data  as  n  -*■“  for  all  t  >0. 

Rapid  surface  vaporization  can  be  modeled  by  means  of 


Knudsen  layer  jump  conditions. 


These  lead  to  the  following 


relationships  between  liquid  surface  properties  (denoted  by  sub¬ 
script  s)  and  vapor  properties  at  the  outer  edge  of  the  Knudsen 
layer  where  equilibrium  is  restored. 

f  =[V1+^?)2  -  *  #  ?]2 

£.  [P  +l)  em  erfc(m)  -  7ij+  \  T  |l  -  me1"2  erfc(m) 

pu/p^  =  1  -  ^2m2  +1-  jl  -  /tt  me™  erfc(m)  ^ 


pu/p, 


-  =  1  -  [ 2m2  +  1 


1  —  ✓'tt  me  er  fc  (m) 


where  y  is  the  ratio  of  specific  heats,  m=u/ /2ET  =  M*'y  /2 ,  and 

12.  Knight,  C. ,  "Theoretical  Modeling  of  Rapid  Surface  Vapori¬ 
zation  with  Back-Pressure,"  AIAA  J.  17,  519  (1979). 
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Figure  2  x-t  Diagram  in  Lab  Reference  Frame 


17 


.  ' AVCO  EVERETT 


M  is  the  flow  Mach  number  at  the  outer  edge  of  the  Knudsen  layer 
at  the  surface.  Thermodynamic  equilibrium  in  the  liquid  will 
be  assumed  so  that  ps  is  plausibly  interpreted  at  the  saturated 
vapor  pressure  corresponding  to  temperature  Tg.  Also  p^ps/R7 
since  the  vapor  will  be  treated  as  a  perfect  gas.  Evaluations 
of  (3)  as  a  function  of  M  are  given  in  Table  1. 

A  Clausius-Cl ayperon  relation  will  be  adopted  to  relate 
saturated  vapor  pressure  and  temperature. 

p*  ■  pb^[^(x  '£;)]  ,4> 


where  p^  =  1  atm,  Tb  is  the  normal  boiling  temperature,  and  Lv 
is  the  (constant)  heat  of  vaporization.  For  aluminum,  y=5/3, 

Tb  =  2730°K,  and  Lv  =  llkJ/g  are  reasonable  choices.  The  Knudsen 
layer  model  embodied  in  (3)  leads  to  a  vapor  stagnation  temper¬ 
ature  which  is  always  less  than  Tg ,  implying  an  energy  flux  into 
the  liquid.  However,  the  temperature  jump  is  small  ( < 1 0 % )  and 
can  be  overlooked.  Thus,  energy  balance  at  the  surface  is  taken 
to  be  expressed  by 


The  absorbed  power  per  unit  surface  area,  dj- ,  is  divided  between 

d 

heat  conduction  into  the  liquid  and  latent  heat  carried  away  by 
the  vapor.  To  close  the  problem,  mass  conservation  requires 

=  pu  (6) 
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K 


M 

P/Ps 

T/Ts 

P/P5 

Pu ,/p  JRT  s 

S  "  2  TT 

0 

1 

1 

1 

0 

0.05 

0.927 

0.980 

0.908 

0.148 

0.1 

0.861 

0.960 

0.827 

0.273 

0.2 

0.748 

0.922 

0.690 

0.465 

0.4 

0.576 

0.851 

0.490 

0.688 

0.6 

0.457 

0.785 

0.358 

0.786 

0.8 

0.371 

0.725 

0.269 

0.817 

1.0 

0.308 

0.669 

0.206 

- -  ,  . 

0.816 

Q  %  MW/cm2 

T  v  °K 
s 

p  'v  a  t  it 

0.1 

2890 

2.07 

1 

3570 

21.8 

10 

4670 

230 

100 

67  50 

2420 

1000 

11,980 

24,400 
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Results  are  best  presented  in  terms  of  nondimensional 
variables.  To  this  end,  let  x  be  the  duration  of  the  pulse 
givinq  rise  to  cr  and  define  £  =  t/x  .  n  =  n/Zot.x  .  fc  = 

a  P  D  S 

xs/xp/ar  f  =  T/Tr,  u  =  u/ar,  p  =  p/pr,  and  p  =  p/pr»  Natural 

choices  for  the  reference  parameters  are  T  =  T,  'v  2730°K, 

r  b 

a  =  / y RtT  ^  1.2  x  105cm/sec,  p_  =  p,  =  1  atm,  and  p  =  p,  /RT, 
r  D  IT  D  FDD 

-  4  3 

'x-  1.2  x  10  g/cm  .  This  leads  to  six  basic  nondimensional 
parameters:  T^T^  y,  L^/RT^ ,  Lv/c£Tb,  Q/a«xp/k£Tb>  and  Prar/ 

p  ^/x^/aC.  where  O  is  a  characteristic  absorbed  power  flux.  Only 
the  case  of  a  top-hat  absorbed  power  flux  has  been  considered 
thus  far — 


q 


a 


O 

0 


0<t<Tp 

t>Tp 


(7) 


Representative  values  of  the  thermal  properties  for  liquid  alumi¬ 
num  are  p  ?  ^  1.6g/cm3,  c^  t  l.lJ/g°K,  t  0.80W/cm°K,  and  ap  'v 

2 

0.45cm  /sec.  Also,  T^  'v  300°K  is  an  interesting  initial  condition. 
These  were  used  in  evaluating  the  nondimensional  parameters  for 
what  follows. 

In  qeneral  the  problem  defined  above  is  coupled  to  the 
exterior  vapor  flow  throuqh  the  Mach  number  at  the  outer  edqe 
of  the  Knudsen  layer,  M.  However,  the  processes  are  decoupled 
so  lonq  as  the  Knudsen  layer  remains  choked  (i.e.,  M  =  1).  This 
special  case  will  be  examined  before  going  on  to  the  fully  coupled 
case.  As  will  be  seen  later,  the  Knudsen  layer  remains  choked 
only  during  a  top-hat  pulse  and  not  afterwards.  Results  for 
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1  usee 


this  special  case  were  generated  numerically  for  xp  = 

and  are  shown  in  Fiqure  3.  It  was  necessary  to  use  an  implicit 

differencing  procedure  for  an  economical  treatment.  Note  that 

the  surface  temperature  approaches  an  asymptote  durinq  the  pulse 
•  2 

for  0  >10  MW/cm  .  That  steady  state  is  defined  by  a  balance  be¬ 
tween  the  vapor  efflux  and  diffusion  terms  on  the  right  side  of 
Eq.  (2).  Note  also  that  this  steady  state  is  reached  increas¬ 
ingly  rapidly  as  O  increases.  The  heat  up  (or  cool  down)  time 

‘  2 

scales  inversely  as  O  ,  as  might  be  expected  from  the  fact  that 
oAp  enters  one  of  the  nondimensional  parameters  above. 

A  straightforward  analysis  leads  to  the  following  steady- 
state  distribution  of  temperature  in  the  liquid. 

T  =  Ti  =  (Tg  -  Ti )  exp  (-xxn/a£)  (8) 

Mass  and  enerqy  balance  at  the  surface  then  leads  to 

<^P_ 

<ct  <TS  -  Tj)  +  Lv]  — i.  =  O  (9) 

/2rRT 

s 

for  a  time  invariant  absorbed  power  O,  where  ^  =  0.816  is  the 
normalized  mass  flux  for  choked  conditions.  (cf.  Table  1) 

This  nonlinear  algebraic  equation  has  been  solved  numerically 
to  qenerate  the  results  in  Table  2.  Note  that  the  saturated 
vapor  pressure  at  steady  state  increases  essentially  linearly 
with  0  whereas  Tg  increases  logarithmically.  The  steady-state 
surface  temperature  is  well  above  critical  conditions 
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ABSORBED  FLUX  (MW/CIVT) 


22 


«  iAVCO  EVERETT 


mr: 


•  2 

(Cf.  Figure  1)  for  0  =  lGW/cm  ,  and  the  predictions  are  highly 
questionable  at  that  flux.  The  model  predictions  are  meaning¬ 
ful  for  aluminum  only  when  the  absorbed  power  flux  q  £  100  MW/cm^. 

Turn  now  to  the  modeling  for  the  vapor  flow  outside  the 
Knudsen  layer.  The  vapor  will  be  treated  as  a  nonabsorbing 
perfect  gas,  as  mentioned  earlier.  It  is  again  convenient  to 
employ  a  moving  reference  plane  attached  to  the  surface.  In 
terms  of  £  =  x  +  x  (t) ,  the  equations  of  motion  become 


If  *  £  (PU)  =  0 


P  U 


0uv  +  = 


H  “  PXs 


pCv  (H  +  U  H)  +  P  If  =  0 


(10) 


In  the  case  of  a  top-hat  energy  input  in  Eq.  (7)  and  a  surround- 
inq  hard  vacuum,  no  shocks  will  appear  anywhere  within  the  flow 
field.  A  method  of  characteristics  then  provides  the  best  ap¬ 
proach  to  solving  for  the  exterior  flow.  Characteristic  equations 


[• 

(10) 

are 

dP 

a 

ds 

df. 

dt 

ii 

dt 

x_  on  -r-r 
s  d  t 

dO 

a 

ds 

dr. 

dt 

yR 

dt 

x  on  i , 
s  dt 

=  u  +  a 


=  u  -  a 


(i  i : 


ds  n  df 
dt  =  0  °n  dt  =  u 
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where  P  =  2a/(y-l)  +  u  and  Q  =  2a/(y-l)  -  u  are  the  Riemann 
invariants  and  s  is  the  entropy.  The  directions  in  the  first 
two  of  these  equations  will  be  referred  to  as  P  and  0  character¬ 
istics,  respectively. 

So  long  as  the  Knudsen  layer  remains  choked  (i.e.,  M  =  1), 
determination  of  surface  conditions  is  independent  of  the  ex¬ 
terior  vapor  flow.  in  that  case  the  external  flow  follows  by 
a  straightforward  marching  procedure  into  £ >  0 ,  starting  with 
known  conditions  at  the  outer  edge  of  the  Knudsen  layer  at  C  =  0 
(The  Knudsen  layer  is  viewed  as  having  infinitesimal  thickness 
in  this  study.)  Note  that  0  characteristics  start  with  zero 
slope  at  £  =  0  and  that  the  slope  will  gradually  become  positive 
in  F,  >0.  The  solution  domain  always  lies  to  the  left  of  the 
terminal  P  characteristic  on  which  a  =  0  and  u  assumes  an  ex¬ 
tremal  value.  The  situation  just  outlined  is  found  to  persist 
so  long  as  the  absorbed  power  flux  is  monotonical ly  increasing 
or  constant  in  time. 

If  it  happens  that  M  <  1,  as  will  happen  when  q  decreases 

the  O  characteristic  will  have  a  negative  slope  and  its  origin 

lies  outside  the  Knudsen  layer.  In  that  case  the  compatibility 

relation  along  the  Q  characteristic  provides  information  deter- 

( 12 ) 

mining  M  and  hence  the  surface  state.  Heat  conduction, 

surface  balance  relations,  and  the  exterior  vapor  flow  are  there 
by  coupled.  Once  the  surface  state  is  established,  the  rest  of 
the  exterior  flow  is  obtained  by  the  same  sort  of  marching  pro¬ 
cedure  noted  in  the  previous  paragraph. 
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c. 


RESULTS  AND  CONCLUSIONS 


Typical  evolution  of  surface  properties  are  displayed  in 

•  2 

Figure  4  for  t  =  1  nsec  and  Q  =  10  MW/cm  .  Taking  an  absorptivity 

a  ~  20%  (as  for  XeF) ,  this  corresponds  to  an  incident  flux  of 

2  2 
50  MW/cm  ,  or  incident  fluence  of  50  J/cm  .  Note  that  the  Mach 


number  at  the  outer  edge  of  the  Knudsen  layer  drops  rapidly  to 
zero  shortly  after  the  pulse-  The  possibility  of  u  <  0  at  the 
surface  is  disallowed  because  condensation  has  not  been  included 


in  the  modeling.  The  surface  temperature  after  the  pulse  is 

significantly  different  from  that  in  Figure  3,  where  M  =  1  at 

all  times  is  assumed.  The  saturated  vapor  pressure  is  also 

shown  below  in  Figure  4.  This  does  not  directly  define  the  force 

acting  on  the  surface.  Neglecting  the  momentum  of  the  liquid, 

2 

(xs>  >  this  is  given  by  the  momentum  flux  of  vapor  from  the 

2  2 

Knudsen  layer:  F=p+pu  =p(l+yM).  The  value  of  F  is 


roughly  55%  of  ps  for  a  choked  Knudsen  layer  and  can  be  evaluated 

using  Table  1  more  generally.  Also  shown  in  Figure  4  is  the 

normalized  mechanical  impulse  delivered  to  the  surface,  Fdt , 

with  F  =  F/p^.  Essentially  all  the  impulse,  when  the  surface 

is  surrounded  by  vacuum,  appears  to  be  delivered  during  the  pulse. 

This  result  will  be  altered  if  plasma  effects  arise,  however. 

Flow  properties  in  the  exterior  flow  at  the  end  of  the 

pulse  (t  =  1)  are  given  in  Figure  5  for  t  =1  psec  and  Q  = 

P 

2 

10  MW/cm  .  Such  information  provides  a  basis  for  assessing 
whether  plasma  formation,  due  to  processes  such  as  laser-induced 
breakdown,  is  likely. 
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Figure  5  Prooerties  in  Exterior  Flow  at  the  End  of  a 

1  nsec  Pulse  for  an  Absorbed  Flux  of  10  MW/cm2 
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An  impulse  coupling  coefficient  C  is  one  parameter 
relevant  to  understanding  damage  mechanisms  associated  with 
pulsed  devices  interacting  with  a  surface.  This  is  normally 


defined  in  terms  of  the  incident  power  density,  ,  but  here 
it  will  be  given  in  terms  of  d[  which  is  typically  a  fraction 
of  .  Thus  the  definition  used  here  is  C 

1  3 

OO  CO 

Ca  =  J  Fdt/Qtp  =  J  Fdt 
o  Q  o 

for  a  top-hat  input.  Values  for  several  conditions  are  shown 

in  Table  3.  The  coupling  coefficient  C  for  the  example  in 

Figures  4  and  5  is  ~7  taps/J  of  absorbed  fluence ,  or~1.4  based 

on  incident  fluence  for  a  ~  20%.  The  value  of  C  is  ~16  at 

a 

2 

at  100  MW/cm  absorbed  and  1  usee,  or  C-3  for  a  ~  20%  (XeF)  and 

2 

C  ~  6  for  ot~40%  (KrF)  at  an  incident  fluence  of  500  J/cm  and 
2 

250  J/cm  ,  respectively. 

In  summary  then,  a  theoretical  model  has  been  constructed 
in  which  two  of  the  primary  mechanisms  influencing  pulsed 
vaporization  of  an  aluminum  surrace  into  surrounding  vacuum 
are  included.  Moreover,  the  results  are  useful  in  assessing  a 
third  possible  process:  absorption  and  transfer  of  energy  in 
the  vapor  phase.  A  summary  of  the  conclusions  follows: 

1)  Generally,  heat  conduction  and  vapor  efflux  represent 
equally  important  mechanisms  for  removing  absorbed 
power  from  the  surface  or  vapor  interface. 
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TABLE  3 


VALUES  OF  IMPULSE  COUPLING  COEFFICIENT  DURING  THE  PULSE 


c 

a  J 

f  1Fdt 
o 

dyne-sec/J 

^ 

lysec 

2ysec 

4ysec 

10  MW/cm2 

7.03 

9.18 

10.81 

30  MW/cm2 

12.25 

13.30 

14.15 

100  MW/cm2 

15.90 

16.90 

— 
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2) 


The  surface  heat-up  time  varies  inversely  with  the 

square  of  the  absorbed  power  and  can  become  small 

compared  to  the  incident  pulse  duration  at  absorbed 
2 

fluxes  >20  MW/cm  ,  or,  for  the  XeF  absorptivity  of 

2 

20%,  at  incident  XeF  fluxes  >100  MW/cm  . 

3)  The  Knudsen  layer  attached  to  the  surface  tends  to  be 
choked  when  the  absorbed  power  is  increasing  and  un¬ 
choked  when  the  power  is  decreasing.  The  problem  is 
fully  coupled  when  the  Knudsen  layer  is  unchoked. 

A  fully  coupled  theoretical  analysis  of  the  sort  pursued  in 
this  Section  is  new  and  has  not  been  published  before  for  the 
case  of  transient  vaporization  into  vacuum. 
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III.  PLASMA  FORMATION  IN  ALUMINUM  VAPOR 

A.  INTRODUCTION 

In  the  previous  section,  an  analysis  was  presented  that 
predicts  aluminum  vapor  conditions  under  irradiation  by  a  high- 
power  UV  laser.  Under  these  conditions,  it  is  possible  for  a 
plasma  to  form.  It  is  the  objective  of  this  section  to  deter¬ 
mine  the  flux  level  for  plasma  onset  at  the  XeF  wavelength 
(353  nm) . 

B .  MODEL 

In  order  to  develop  a  tractable  model  for  plasma  formation 
in  aluminum  vapor  irradiated  by  an  XeF  laser,  simplifying  assump¬ 
tions  must  be  made  about  the  electronic  structure  of  aluminum. 

We  use  the  energy  level  diagram  of  aluminum  shown  in  Figure  6 
as  the  starting  point  for  our  approximations.  This  figure  shows 
all  of  the  excited  states  up  to  4.7  eV  above  the  ground  state  and 
the  first  ionization  limit  at  5.98  eV.  The  first  excited  state 
is  3.14  eV  above  the  ground  state  and  is  optically  accessible 
from  the  ground  state  and,  thus,  is  expected  to  be  important  in 

energy  storage  in  aluminum  vapor.  The  second  excited  state  is 

2  4 

a  (3s3p)  P  state.  Since  optical  transition  to  this  state  from 
the  ground  state  is  both  spin  and  angular  momentum  forbidden, 
its  influence  on  energy  storage  in  the  vapor  is  probably  negli¬ 
gible  and,  thus,  has  been  omitted  from  our  model.  Immediately 
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Figure  6  Energy  Level  Diagram  for  Aluminum 


above  the  state  are  the  (3s^3d)^D  and  (3s^4p)^P  states  which 

lie  within  kT  of  each  other  for  all  temperatures  of  interest  for 

2 

our  application.  The  next  excited  state,  a  S  state,  occurs  at 

4.67  eV  and  at  higher  energies,  excited  states  are  plentiful 

up  to  the  first  ionization  limit. 

For  our  calculations,  we  have  assumed  that  the  electronic 

manifold  of  staes  in  aluminum  can  be  represented  by  the  ground 

state,  Al,  two  excited  states,  Al*  and  Al**,  and  a  singly  ionized 

state,  Al+.  One  excited  state,  Al*,  corresponds  to  the  first 

2  2 

excited  state,  (3s  4s)  S  of  aluminum.  Excitation  to  this  state 

is  expected  to  be  the  rate  limiting  step  for  ionization,  since 

this  process  requires  the  transfer  of  3.1  eV  of  energy  and  only 

2.9  eV  are  required  to  complete  the  ionization.  The  other  ex- 

2  2 

cited  state,  Al**,  represents  a  combination  of  the  (3s  3d)  P 
states.  This  state  can  be  easily  formed  by  transitions  from 
both  the  ground  and  first  excited  states.  Once  the  Al**  state 
is  formed,  excitation,  and  subsequent  ionization,  is  expected 
to  occur  rapidly  due  to  the  density  of  states  between  the  Al** 
state  and  the  ionization  limit,  therefore,  including  more  highly 
excited  states  in  our  model  is  not  considered  necessary.  Only 
the  singly  charged  ion  is  included  in  our  model  since  the  mul¬ 
tiply  ionized  states  do  not  become  important  until  temperatures 
of  a  few  eV  are  reached. 

In  our  model,  excitation  of  aluminum  vapor  is  assumed  to 
occur  by  two  mechanisms;  electron  impact  excitation  and  photo- 
ionization.  These  processes  can  be  represented  by  the  following 
reactions : 
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Electron  impact  excitation  and  ionization: 


Al  +  e  -+  Al*  +  e  (1) 

A1  +  e  ->  Al  *  *  +  e  (2) 

Al*  +  e  -*•  Al**  +  e  (3) 

Al**  +  e  Al+  +  2e  (4) 

Photo ionization 

Al*  +  hv  -*■  A1+  +  e  (5) 

Al**  +  hv  -*■  Al+  +  e  (6) 


Direct  electron  impact  ionization  of  the  Al  and  Al*  states  is  not 
included  in  our  model  since  it  would  be  much  slower  than  the 
competing  process  of  excitation  to  Al**  and  subsequent  ioniza¬ 
tion  from  that  state. 

Other  processes  which  may  affect  the  electron  temperature 
and,  hence,  the  rate  of  vapor  heating  are  momentum  transfer  from 
the  electrons  to  the  heavy  particles  (reaction  (7))  and  inverse 
Bremsstrahlung  absorption  of  the  laser  radiation  by  electrons 
in  the  presence  of  neutral  or  ionic  aluminum. 

Also  necessary  for  our  modeling  of  plasma  formation  in 
aluminum  vapor  are  rate  constants  for  the  reactions  listed  above. 

A  first  cut  at  the  electron  impact  excitation  rate  constants  for 
reactions  ( 1 )  —  ( 3 )  can  be  made  by  using  the  following  expression:11 

f  -W  /kT  /W  \ 

k  .  =  1.58  x  lO-3  - - —  e  e  e  Pi  )  (cm3/ sec) 

13  sr  w  V  V 

e  e 

where  f  is  the  optical  oscillator  strength  of  the  transition  i»j, 
We  is  the  energy  of  the  transition,  is  the  electron  temperature 
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and  P (x)  is  a  tabulated  function. 


(13) 


The  rate  constant  for  the 


reverse  reaction  can  be  obtained  by  detailed  balancing. 


W  /kT 
e  e 


where  g^  and  g^  are  the  degeneracies  of  the  lower  and  upper  states, 
respectively . 

A  reasonable  rate  constant  for  the  remaining  electron  impact 
process,  Al**  +  e  ■*  Al+  +  2e,  (reaction  (4))  is  more  difficult  to 
determine.  While  a  rate  constant  for  direct  ionization  from  the 
Al**  state  can  be  calculated,  this  constant  will  be  a  lower 

bound  on  the  true  value  since  ionization  of  the  Al**  state  can 
also  occur  via  multistep  excitation  followed  by  ionization  from 
a  highly  excited  electronic  state  whose  rate  constant  for  ioniza¬ 
tion  will  be  much  larger  than  that  for  the  Al**  state.  This 
multistep  process  will  increase  the  effective  rate  for  ionization 
and,  indeed,  will  probably  be  the  dominant  ionization  channel. 

In  our  model,  we  have  therefore  assumed  that  the  lifetime  for 
ionization  from  the  Al**  state  is  approximately  three  times 
longer  than  the  lifetime  for  excitation  to  the  next  more  ener¬ 
getic  level,  since  three  transtions  of  this  energy  separation 
are  just  sufficient  for  ionization.  The  rate  constant  obtained 
using  this  approximation  is  appreciably  larger  than  that  for 
direct  ionization. 

13.  Allen,  C.W.,  Astrophysica 1  Quantities,  2nd  Edition,  Athlone 
Press,  London,  1963,  p.  42. 
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As  will  be  seen  from  our  computer  results,  the  time  nec¬ 
essary  for  plasma  formation  is  sensitive  to  the  rate  constant 
for  the  reverse  of  reaction  (4) .  Therefore,  a  more  accurate 
method  of  estimating  this  rate  constant  would  be  desirable  and 
would  increase  the  reliability  of  our  computer  predictions. 

The  expression  given  above  for  the  electron  impact  excita¬ 
tion  rate  constants  was  derived  using  a  Maxwellian  energy  distri¬ 
bution  for  the  electrons.  For  appreciable  excited  state  densities 
which  can  occur  in  the  regime  of  interest  for  this  study,  this 
approximation  breaks  down  and  a  more  accurate  method  of  calculat¬ 
ing  rate  constants  which  allows  for  a  non-Maxwe Ilian  distribution 
should  be  considered  in  future  work. 

2  2 

Cross  sections  for  photoionization  from  the  4s  S,  3d  D  and 

2 

4p  P  states  which  appear  in  our  model  (reactions  (5)  and  (6) )  can 

be  calculated  using  the  quantum  defect  method  of  Burgess  and 
( 14 ) 

Seaton.  The  cross  section  for  ionization  from  the  4s  state 

possesses  a  Cooper  minimum.  This  phenomenon  has  also  be  observed 
for  ionization  from  s  states  in  the  alkalis^5^  and  rare  gases. 

The  cross  section  for  ionization  by  an  XeF  photon  will  occur  near 
this  Cooper  minimum,  thus  causing  that  process  to  have  an  ano- 
nalously  small  cross  section.  No  such  minimum  occurs  in  the  3d 
cross  section  and  while  a  minimum  does  occur  in  the  4p-*-3d  cross 

14.  Burgess,  A.  and  Seaton,  M.J.,  Mon.  Notices  Roy.  Astron. 

Soc.  120,  121  (1960). 

15.  Fano,  U. ,  Cooper,  J.W.,  Rev.  Mod.  Phys.  £0,  441  (1968). 

16.  Duzy,  C.  and  Hyman,  H.A.,  Phys.  Rev.  A.  22,  1878,  (1980). 
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section,  ionization  by  an  XeF  photon  is  not  near  this  dip  in  the 
cross  section.  For  absorption  of  353  nm  photons  the  cross  sec¬ 
tions  for  the  4s,  3d,  4p  states  are  4.9  x  10-19  cm2,  1.1  x  10-17  cm 
—  1 8  2 

and  7.0  x  10  cm  ,  respectively.  In  our  model,  the  Al*  state  has 

- 1  9  2 

a  photoionization  cross  section  of  4.9  x  10  cm  while  that  for 

_  1 8  2 

the  Al**  state  is  9.4  x  10  cm  .  This  difference  in  behavior 

of  the  4s  state  toward  the  incident  radiation  is  the  main  reason 

for  our  choosing  to  include  two  excited  states  in  our  model. 

( 5 ) 

This  approach  differs  from  that  used  by  Rosen,  et  al. 

Of  the  remaining  processes  mentioned  previously  and  included 
in  our  model,  neutral  inverse  Bremsstrahlung  (reaction  (9)  in 
our  code)  is  expected  to  have  the  greatest  effect  on  plasma  form¬ 
ation.  All  calculations  of  the  cross  section  for  this  process 
to  date  have  been  based  on  the  assumption  that  absorption  of  a 

photon  by  an  electron  causes  only  a  negligible  change  in  the 

(7) 

momentum  of  the  electron.  This  assumption  is  obviously  not 

valid  when  absorption  of  a  3.5  eV  XeF  photon  is  being  considered. 

However,  until  a  more  appropriate  estimate  is  available,  we  will 

..  .  .  .  (17) 

use  the  existing  expression 

a  =  1.88  x  10~15  /a  A3T  (cm5) 

e 

where  a  is  the  polarizability  of  the  atom  in  cm3,  \  is  the  wave¬ 
length  of  the  absorbed  radiation  in  cm  and  Tg  is  the  mean  electron 
temperature  in  eV. 


17.  Hyman,  H.A.,  Kivel,  B.  and  Bethe,  H.A.,  AERL  Report  #AMP 
377,  April  1973. 
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Inverse  Bremsstrahlung  absorption  in  the  presence  of  an 

(18) 

atomic  ion  (reaction  (8))  has  the  following  cross  section: 

o  =  1.25  x  10"25  X3//?^  (cm5) 

where  X  and  T  are  as  defined  above, 
e 

The  reactions  described  above  and  their  rate  constants  or 
cross  sections  were  incorporated  into  a  computer  code  and  the 
time  evolution  of  the  system  was  followed.  The  starting  densi¬ 
ties  of  the  various  species  were  obtained  from  the  calculation 
of  vapor  pressure  and  temperature  described  in  Section  II.  Since 
the  plasma  code  assumes  constant  vapor  density,  we  chose  as  ini¬ 
tial  conditions  values  just  outside  the  Knudsen  layer,  at  a  time 
of  1/2  the  pulse  length.  As  will  be  shown  below,  the  choice  of 
this  particular  time  within  the  pulse  is  not  important  because 
plasma  breakdown  occurs  at  fluxes  high  enough  that  the  vapor 
reaches  steady  state  early  in  the  pulse. 

C .  RESULTS 

We  will  now  present  numerical  results  from  the  computer 
model  described  avove .  The  initial  conditions  for  our  computer 
model  were  obtained  as  described  previously  using  the  results  of 
Section  II.  In  Table  4,  the  heavy  particle  density  and  the  heavy 
particle  and  electron  temperatures  are  listed  as  a  function  of 
incident  flux.  The  remaining  densities  can  be  generated  usinn 
this  information  and  the  thermal  equilibrium  assumption.  In  the 


Zel'dovich,  Ya  B.,  Raiser,  Yu  P. ,  Physics  of  Shock  Waves 
and  High-Temperature  Hydrodynamics  Phenomena,  Vol .  1, 
p.  259,  Academic  Press,  New  York,  1966. 
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previous  section,  the  results  were  expressed  in  terms  of  the 
absorbed  flux,  since  that  is  the  dominant  parameter  in  that  cal¬ 
culation.  In  this  section  the  results  are  expressed  in  terms 
of  the  incident  flux,  which  is  the  more  important  parameter 
here.  The  two  are  related  by  the  absorptivity  a,  which  is  taken 

to  be  20%  at  the  XeF  wavelength.  Beginning  at  the  highest 

2  2 

incident  flux  density  considered  (250  MW/cm  )  (or  50  MW/cm 
absorbed) ,  we  see  in  Figure  7  that  the  electron  temperature 
rises  to  1.0  eV  in  0.057  ysec,  at  which  time  breakdown  may  be 
considered  to  have  occurred  due  to  the  rapidity  with  which  ions 
are  formed  at  this  temperature  and  the  particle  densities  con¬ 
sidered  in  this  work.  Above  1  eV  our  model  with  its  assumed 
Maxwellian  electron  distribution,  becomes  increasingly  unreal¬ 
istic,  so  the  calculation  is  stopped  at  Tg  =  1.5  eV. 

Figure  7  shows  that  the  temperature  of  the  heavy  particles, 

—  8 

Tg,  also  begins  to  rise  rapidly  after  10  sec,  and  momentum 
transfer  is  sufficient  to  equilibrate  energy  between  electrons 
and  atoms. 

The  electron  density  (n  )  shows  a  similar  history,  as  seen 
in  Figure  8.  The  concentration  of  excited  aluminum  atoms  also 
rises  rapidly  during  the  breakdown  period. 

2 

At  an  incident  flux  density  of  100  MW/cm  the  breakdown 

2 

occurs  more  slowly  than  at  250  MW/cm  .  As  we  see  in  Figure  9, 

the  electron  temperature  reaches  1.0  eV  in  0.092  sec.  Figure  10 

shows  an  avalanche  of  ionization  and  atomic  excitation,  also 

2 

delayed  relative  to  the  250  MW/cm  case. 

40 


^1AV CO  EVERETT 


TEMPERRTURE  (EV) 


J8047 


""IjqM . . 


oJ,/  ""'ll 

TIME 


(SEC) 


TTTT 


1j  0t»  r  r  ",tTj  qi*t,t""'1  01? " '  ""'i  o'* 


Figure  7  Temperature  of  Electrons  (Te)  and  Temperature  of 
Aluminum  Atoms  and  Ions  (T„)  for  an  Incident  Flux 
Density  of  250  MW/cm^ 
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Concentration  of  Electrons  (ne)  and  Excited  Aluminum 
Atoms  ([A2,*]  and  [A£**])  at  250  MW/cm^ 
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Figure  10 


Electron  and  Excited  Atom  Densities  at  100  MW/cm 


2 
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Rates  of  individual  reactions  (1)— (6) (8)  and  (9)  are  plotted 

3 

in  Figures  11  and  12  in  part/cm  -sec  for  an  incident  flux  of 
2 

100  MW/cm  .  Reactions  going  in  the  forward  direction  appear  in 
Figure  11,  wh.i  le  those  going  backwards  are  shown  in  Figure  12. 

We  observe  in  these  figures  three  distinct  reactions  regimes, 
(a)  During  the  first  10  sec  the  production  of  excited  aluminum 
atoms  dominates  via  the  reactions 


Al 

+ 

e 

-*■  Al*  +  e , 

(1) 

Al 

+ 

e 

-*■  Al**  +  e. 

(2) 

Al** 

+ 

e 

Al+  +  2e. 

(4-) 

In  this  regime  the  ion  density  falls  slightly.  This  is  associated 

with  the  slight  decrease  in  electron  temperature  which  occurs  as 

-10  -9 

the  aluminum  is  excited.  (b)  From  2  x  10  to  4  x  10  the  scene 
is  relatively  quiet.  The  chief  reaction  is 

Al**  +  hv  -*•  Al+  +  e,  (6) 

but  almost  all  of  the  ionization  is  neutralized  by  reaction  (4-). 

Reaction  (6)  gradually  builds  up  the  electron  concentration  so 

that  reaction  (2)  comes  to  life  again.  In  addition  to  creating 

more  electrons,  reaction  (6)  is  also  the  main  mechanism  for 

heating  electrons.  Since  the  Al**  state  is  only  1 . 9  eV  below 

the  ionization  limit,  absorption  of  an  XeF  (3.5  eV)  photon 

creates  an  electron  with  1.6  eV  of  energy,  which  is  considerably 

above  the  mean  electron  energy.  The  importance  of  this  will  be 

-  9 

evident  later.  (c)  After  7  x  10  sec,  the  combination  of  re¬ 
actions  (2)  and  (6)  leads  to  avalanche  ionization. 
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Reverse  Reaction  Rates  at 


100  MW/cm2 
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2 

The  breakdown  mechanism  is  similar  at  50  and  250  MW/cm  . 

2  2 

At  an  incident  flux  of  50  MW/cm  (10  MW/cm  absorbed)  we  find 
that  breakdown  is  sensitive  to  the  rate  constant  for  the  reverse 
of  reaction  (4).  The  above  discussion  of  our  model  indicated 
uncertainties  concerning  this  reaction.  The  recombination  rate 
constant  (4-)  is  obtained  by  multiplying  the  ionization  rate 
constant  by  the  Saha  expression  in  which  the  energy  W  of  the 
reaction  Al**  +  e  -*■  A+  +  2e  (4)  appears  in  a  positive  exponent. 
Such  a  formalism  implies  a  one-step  process.  Actually,  recom¬ 
bination  will  involve  several  steps.  To  model  this,  we  replace 
the  energy  in  the  exponent  by  W/B,  with  B  an  adjustable  parameter. 
Thus,  8=1  corresponds  to  a  single  step  for  reaction  (4-1), 
while  6=3  corresponds  to  three  steps. 

The  cases  presented  above  are  for  B  =  1.  At  50  MW/cm2 
incident,  we  have  breakdown  (Tg  -*•  1.0  eV)  in  0.21  usee  if  B  =  1, 
as  displayed  in  Figures  13  and  14.  If  6  =  2  or  3,  however,  no 
breakdown  occurs  within  2  ysec  (Figures  15  and  16) . 

Since  reaction  (4-)  is  the  only  mechanism  for  removal  of 
electrons  from  the  vapor,  one  would  expect  a  decrease  in  its 
rate  constant  (i.e.,  increasing  8)  to  decrease  the  time  necessary 
for  breakdown.  Actually,  the  opposite  occurs. 

This  seemingly  paradoxical  situation  can  be  explained  as 
follows.  The  explanation  hinges  on  the  fact  that  the  main  elec¬ 
tron  heating  mechanism  in  our  model  is  photoionization  from  the 
Al**  state.  As  the  rate  constant  for  recombination  is  increased 
(6  decreased)  more  Al**  states  are  formed  relative  to  the  number 
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Figure  13  Electron  and  Gas  Temperatures  at  50  MW/cnT 
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15  Electron  and  Gas  Temperatures  at  50  MW/cm  . 

Here  3  =  2;  i.e.,  the  rate  constant  for  reaction 
(4-)  is  reduced.  Breakdown  occurs  after  3.1  psec. 
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Electron  and  Excited  Atoms  at  50  MW/cm 
Here  8  =  2. 
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of  ions  present.  These  Al**  states  can  photoionize  which  pro¬ 
duces  hot  electrons  which  in  turn  increase  the  mean  electron 


temperature.  This  higher  electron  temperature  increases  the 
rate  of  excitated  state  production  via  electron  impact,  so  more 
Al**  states  are  formed  and  the  process  continues.  If  the  num¬ 
ber  of  Al**  states  available  for  photoionization  is  decreased 
by  decreasing  the  recombination  rate,  the  electron  heating  rate 
will  decrease  and  the  onset  of  breakdown  will  be  delayed  as  is 
seen  in  our  results.  in  other  words,  to  achieve  breakdown,  it 
pays  to  sacrifice  a  few  electrons  in  order  to  heat  up  the  re¬ 
mainder. 

The  effect  of  6  on  the  time  to  achieve  Tg  =  1.0  eV  is 
summarized  in  Table  5.  it  appears  that  the  rate  constant  for 
reaction  (4-)  has  a  critical  effect  on  whether  or  not  breakdown 
occurs  within  a  given  time. 

These  results  naturally  depend  upon  the  initial  conditions. 
The  time  to  achieve  breakdown  evidently  decreases  as  the  initial 
electron  density  increases. 

2 

No  breakdown  at  all  occurred  at  25  MW/cm  incident. 

D.  CONCLUSIONS 

2 

Plasma  formation  was  not  predicted  at  25  MW/cm  for  any  of 

the  sets  of  rate  constants  considered  in  this  work  while  at 
2 

100  MW/cm  plasma  was  formed  very  rapidly  regardless  of  the 
assumptions  used  for  determining  the  rate  constants.  For  flux 
levels  between  these  two  extremes  our  results  indicate  that  the 
relative  magnitude  of  the  photoionization  and  electron  impact 
ionization  rates  determine  whether  a  plasma  is  present. 
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(4  5) 

For  a  0.5  usee  pulse  for  which  data  exist  '  our  code 

results  suggest  that  a  plasma  should  not  be  present  for  fluences 

2  2 
<  10  J/cm  and  will  almost  certainly  be  present  above  5Q  J/cm  . 

This  conclusion  is  consistent  with  the  existing  experimental 

results.  Rosen,  et  al.^  first  observed  a  decrease  in  reflected 

2 

laser  intensity  at  ~30  J/cm  .  This  value  lies  within  the  ex- 
prected  range  for  the  onset  of  plasma  formation. 

Several  improvements  could  be  made  in  our  model  which 
would  increase  its  reliability  and  we  suggest  their  inclusion 
in  future  work.  The  inclusion  of  a  non-Maxwell ian  electron 
distribution  for  determining  the  electron  impact  rate  constants 
is  one  such  improvement.  This  correction  could  be  very  important 
when  large  numbers  of  excited  atoms  are  present.  Further  thought 
should  also  be  given  to  the  expression  used  for  estimating  elec¬ 
tron  impact  ionization. 

Uncertainty  also  exists  in  our  calculation  of  the  cross 
section  for  photoionization  from  the  4p  state  due  to  the  diffi¬ 
culty  in  estimating  the  quantum  defect  for  the  4p-*-3d  transition. 

A  more  detailed  calculation  could  eliminate  this  ambiquity. 

One  effect  which  was  not  considered  in  our  model  is  re¬ 
radiation  by  the  aluminum  vapor  plasma.  While  the  net  energy 
lost  by  this  mechanism  is  probably  small  when  averaged  over  the 
entire  vapor,  reradiation  is  the  mechanism  by  which  energy  is 
coupled  from  the  plasma  to  the  surface,  and  therefore  should  be 
studied  in  detail.  Multiply  -  charged  ions  should  also  be  included. 
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APPENDIX  A 


UV/VISIBLE  LASER  EFFECTS 

1 .  INTRODUCTION 

While  a  great  deal  of  work  has  been  done  on  laser  effects 

with  pulsed  infrared  lasers, ^  only  recently  have  UV  lasers  of 

12) 

suitable  energy  levels  been  developed.  In  the  IR,  the  inter¬ 

action  between  a  pulsed  laser  and  an  aluminum  surface  in  air  is 
dominated  by  the  formation  of  a  plasma,  which  totally  absorbs 
the  laser  energy.  The  plasma  then  reradiates  a  large  fraction  of 
this  energy  in  the  UV.  Since  aluminum  is  a  very  good  reflector 
in  the  IR,  but  a  good  absorber  in  the  UV,  this  leads  to  higher 
surface  coupling  than  would  be  expected  if  the  plasma  were  absent. 

In  addition,  the  plasma  generates  a  significant  impulse  level  on 
the  target  by  remaining  at  high  pressure  for  a  significant  period 
of  time.  With  a  UV  laser,  where  the  photons  are  much  more  ener¬ 
getic  and  where  the  basic  absorptivity  of  the  surface  is  higher, 
the  phenomenology  may  be  different. 

We  measured  the  thermal  coupling  of  the  pulsed  KrF/XeF  laser 
beam  to  the  target  and  the  impulse  generated  by  the  laser  using 

standard  techniques^  (thermocouples  and  a  pendulum,  respectively). 

( 2 ) 

The  laser  and  its  diagnostics  have  been  described  previously. 

For  KrF  experiments  it  was  operated  so  as  to  produce  pulses  of 

500  nsec  pulse  length  and  up  to  20  J  energy.  A  stable  resonator 

T~.  McKay,  J.A.,  Schriempf,  J.T.,  Cronburg,  T.L.,  Eninqer,  J.E., 
and  Woodroffe,  J.A.,  Appl.  Phys.  Lett.  3^,  125  (1980). 

2.  Jacob,  J.  H. ,  Hsia,  J.  C.,  Mangano,  J.  A.,  and  Rokni ,  M. , 

J.  Appl.  Phys.  50,  5130  (1979). 
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was  used,  leading  to  multi-mode  output.  The  fluence  level  on 
the  targets  was  varied  by  locating  targets  at  different  positions, 
relative  to  a  50  cm  focal  length  lens.  The  environment  was  1  atm 
air.  The  minimum  spot  diameter  used  was  0.5  cm.  The  targets 
were  made  of  2024-T3  aluminum,  0.05  cm  thick  (0.020  in.).  The 
pendulum  weight  was  7.5  g. 

For  most  of  the  XeF  work,  a  laser  pulse  length  of  600  nsec 
was  used,  although  some  of  the  impulse  coupling  data  were  obtained 
using  pulse  lengths  up  to  1.5  usee.  Output  energies  up  to  7  J 
were  used.  The  XeF  experiments  were  performed  in  a  vacuum 
chamber  which  was  placed  at  the  focal  point  of  the  50  cm  focal 
length  lens  which  focussed  the  laser  output.  By  moving  the  lens, 

2 

the  spot  size  on  our  targets  was  varied  between  ~0.04  and  ~0.6  cm  . 
We  used  a  pendulum  with  a  15  cm  arm  and  a  mass  of  ~1.2  gm. 

2 .  RESULTS 

Figure  A-l  shows  the  KrF  results  for  thermal  coupling.  We 
have  plotted  the  fraction  of  the  incident  laser  energy  that  remains 
in  the  target  after  the  pulse  is  over  versus  incident  fluence.  The 
advantage  of  this  procedure  is  that  the  energy  out  of  the  laser  is 
known  more  accurately  than  the  fluence  distribution  over  the  spot. 
As  long  as  the  spot  is  not  too  small,  the  difference  will  be  small 
between  coupling  measured  in  this  way,  and  that  measured  by  re¬ 
lating  the  temperature  rise  in  the  center  of  the  spot  to  the 
fluence  incident  on  the  center  of  the  spot.  The  fraction  of 
energy  coupled  is  in  the  range  40%.  The  low  flux  room  temperature 
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FRACTION  OF  INCIDENT  PULSE  ENERGY 
DEPOSITED  ON  TARGET 


jziee  INCIDENT  FLUENCE  (  J/CM2  ) 


Figure  A-l  Thermal  Coupling  of  KrF  Laser  to  Bate  Aluminum. 

Pressure:  1  atm  pulse  energy:  7-10  J  pulse  length: 

0.5  usee.  The  heavy  line  is  a  guide  to  the  eye 
through  the  data.  The  light  lines  are  analytical 
bounds  on  the  coupling,  based  on  40%  inherent  ab¬ 
sorptivity;  a  is  the  thermal  diffusivity  of  aluminum. 
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The 


absorptivity  at  this  wavelength  is  typically  30-40%. ^ 

phenomenological  behavior  is  shown  at  the  bottom  of  the  plot; 

2 

melt  is  visible  on  the  targets  at  fluences  >6  J/cm  ;  at  fluences 
~  2 

>20  J/cm  ,  a  crater  rim  is  observed,  indicating  that  melt  has 
been  splashed  to  the  edge  of  the  spot,  but  no  mass  loss  from  the 
target  was  observed  (to  a  level  of  0.2  mg)  at  any  fluence.  We 
have  also  shown  in  Figure  A-l  two  theoretical  bounds  on  the 
coupling,  both  related  to  energy  removal  by  metal  vapor  flow 
from  the  target  surface.  A  lower  bound  can  be  obtained  by  as¬ 
suming  that  as  soon  as  the  1  atm  vapor  pressure  is  achieved,  all 
of  the  extra  incident  energy  is  carried  off  by  the  vapor.  An 
upper  bound  is  obtained  by  assuming  that  no  energy  is  carried 
off  until  steady-state  vaporization  is  achieved  (at  a  much 
higher  vapor  pressure  aned  surface  temperature) .  The  value  for 
steady-state  temperature  was  obtained  using  the  model  of  Ref.  4. 

The  data  can  be  seen  to  fall  off  between  the  two  bounds. 

Figure  A-2  contains  the  results  of  the  thermal  coupling 
measurements  with  XeF.  We  have  plotted  the  fraction  of  the 
incident  laser  energy  that  remains  in  the  target  after  the  pulse 
is  over  versus  the  incident  fluence.  The  thermal  coupling  at 
moderate  fluences  is~20%.  The  low  fluence  limit  of  ~22%  in¬ 
dicated  in  this  figure  was  obtained  from  one  of  our  samples 
after  polishing  using  a  Cary  17  Spectrophotometer.  This  is  in 

3.  Touloukian,  Y.  S.  and  Ho,  C.  Y.,  eds.,  Thermophysical  Properties 
of  Selected  Aerospace  Materials,  (TEPIAC7CINDAS ,  Purdue  Univer¬ 
sity)  (1976). 

4.  Knight,  C.,  A1AA  J.  17,  519  (1979). 
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FRACTION  OF  INCIDENT  ENERGY 
DEPOSITED  IN  TARGET 


LOW  FLUX 
ABSORPTIVITY 


.i  . 

iV  • 


•  PRESSURE  =  lOO/i 
O  PRESSURE  *  I  ATM 


|- MELT 


(-♦MELT  SPLASH 


INCIDENT  FLUENCE  (J/cm ) 


Figure  A-2  Thermal  Coupling  of  XeF  to  an  Aluminum  Surface  as 
a  Function  of  Incident  Fluence.  A  pulse  length  of 
0.6  s  and  energies  up  to  7  J  were  used.  Results  in 
air  at  0.1  torr  (dots)  and  760  torr  (open  circles) 
are  shown.  The  low  fluence  absorptivity  was  obtained 
f rom  Ref.  3 . 
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agreement  with  previously  reported  values  for  2024  aluminum  of 
(3) 

20-30%.  While  most  of  our  data  were  obtained  at  a  background 

pressure  of  0.1  torr,  three  shots  were  taken  at  atmospheric  pres¬ 
sure  and  are  indicated  in  Figure  A-2.  Both  in  air  and  vacuum, 

2 

at  fluences,^  4  J/cm  ,  melt  was  observed  on  the  surface  of  our 

2 

targets.  For  f luences Z  100  J/cm  a  "crater  rim"  caused  by  melting 
aluminum  being  splashed  to  the  edges  of  the  laser  spot  was  ob¬ 
served.  It  should  be  noted  that  no  melt  splash  was  present  for 

2 

fluences  up  to  30  J/cm  ,  even  though  melt  splash  was  observed 
.  2 

with  KrF  in  air  at  fluences  >20  J/cm  .  This  is  consistent  with 

the  measured  thermal  coupling. 

All  of  our  measurements  were  obtained  at  fluence  levels 
2  2 

(£4  J/cm  )  and  flux  levels  (^7  MW/cin  )  sufficient  to  cause  sur¬ 
face  melt,  thereby  causing  our  samples  to  behave  like  polished 
surfaces.  This  is  consistent  with  the  fact  that  our  measured 

absorptivities  are  similar  to  the  low  fluence  absorptivity  of  a 

(3) 

polished  surface.  However,  for  very  small  incident  fluences 

where  no  surface  melting  occurs,  the  roughness  of  the  target 

surface  may  cause  an  increase  in  the  observed  absorptivity. 

Figure  A-3  shows  the  results  for  impulse  generated  by  the 

KrF  laser  versus  incident  fluence.  We  have  plotted  impulse  on 

the  left-hand  side  and  a  coupling  coefficient  C  for  impulse 

2 

versus  energy  on  the  right-hand  side.  At  0.3  J/cm  ,  no  surface 
flash  or  pendulum  motion  were  observed.  For  all  the  other  points, 
a  surface  flash  was  seen  (even  though  at  the  lower  fluences 
surface  melt  was  not  apparent).  In  the  infrared,  it  is  well 
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Figure  A-3  Impulse  Produced  by  a  KrF  Laser  on  Aluminum.  Pres¬ 
sure:  1  atm  pulse  energy:  20  J,  pulse  length:  0.5 
usee.  The  heavy  line  is  a  guide  to  the  eye  through 
the  data.  The  light  lines  are  calculations  for 
plasma  (laser  supported  detonation  wave)  impulse 
with  an  ignition  threshold  of  3  J/cm2  and  vapor 
induced  impulse  with  40%  inherent  absorptivity. 
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known  that  a  surface  plasma  can  form  which  can  lead  to  high 

thermal  coupling  and  to  significant  pressure  and  impulse  levels. 

We  have  taken  a  model  developed  for  plasma  impulse  coupling  by 

an  IR  laser  via  a  laser  supported  detonation  (LSD)  wave ^  and 

used  it  without  modification  for  the  fluxes  and  pulse  length  of 

these  experiments,  assuming  a  plasma  ignition  threshold  level  of 
2 

3  J/cm  .  The  calculation  is  shown  on  Figure  A-3.  For  comparison 

we  have  also  done  a  calculation  for  the  impulse  caused  by  steady- 

state  bulk  vaporization  recoil.  The  timescale  for  steady- 

state  vaporization  to  be  reached  is  longer  than  the  pulse  length 

for  the  entire  range  of  experimental  fluences.  The  vaporization 

calculation  is  based  on  a  steady-state  condition  being  reached, 

and  is  thus  an  upper  bound.  This  becomes  a  more  accurate  descrip- 

2  2 

tion  as  the  fluence  increases  towards  10  J/cm  .  The  general 
levels  of  the  two  calculations  are  in  the  right  range,  though 
neither  is  accurate,  suggesting  that  a  more  detailed  model  of  the 
physics  will  be  required. 

Figure  A-4  shows  the  results  of  our  impulse  coupling 
measurements  for  XeF.  The  impulse  coupling  coefficient  (dyne- 
sec/joules)  is  plotted  versus  incident  fluence.  Results  for 
background  pressures  of  0.1  and  2  torr  are  shown;  there  is  no 
discernable  difference  between  the  two  sets  of  results.  Measure¬ 
ments  to  the  impulse  coupling  coefficient  were  also  performed  at 

5^  Reilly,  J.  P.,  Ballantyne,  A.,  Woodroffe,  J.  A.,  AIAA.  J.  17, 
1098  (1979). 

6.  Shui,  V.  H.,  Young,  L.  A.,  Reilly,  J.  P.,  AIAA  J.  16,  649 
(1978) . 
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Figure  A-4 


Impulse  Coupling  Coefficient  for  XeF  to  an  Aluminum 
Surface  as  a  Function  of  Incident  Fluence.  Energies 
of  0. 8-2.0  J  were  used. 
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atmospheric  pressure  and  75  torr  and  these  results  are  given 
in  Figure  A-5.  Again,  the  results  at  1  atm  and  75  torr  exhibit 
similar  behavior.  It  should  be  noted  that  no  sudden  disappearance 
of  impulse  coupling  as  a  function  of  decreasing  fluence  was  ob¬ 
served.  Instead,  the  pendulum  motion  gradually  decreased  as  the 

2 

fluence  decreased  until  at  ~4  J/cm  ,  any  motion  became  too  small 
for  us  to  detect.  At  all  fluences  where  pendulum  motion  was  ob¬ 
served,  a  surface  flash  was  also  present  and,  moreover,  this  sur¬ 
face  flash  persisted  until  the  fluence  was  decreased  below  ~2  J/cm 

We  observed  this  behavior  even  though  the  target  surface  had 
been  melted  by  a  number  of  previous  shots.  We  believe  it  is  due 
to  the  presence  of  surface  defects.  Scans  of  the  surface  by  SEM, 
in  fact,  show  a  large  number  of  defects  on  pre-melted  surfaces, 
as  well  as  on  "as-received"  surfaces. 

3.  CONCLUSION 

In  summary,  we  have  observed  high  thermal  coupling  of  KrF 
and  XeF  laser  radiation  to  an  aluminum  surface  (~40%  and  ~20%, 
respectively) ,  and  the  appearance  of  significant  impulse  levels 
at  low  fluence. 
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IMPULSE  COUPLING  COEFFICIENT  (DYNE-SEC/J) 


Figure  A-5  Impulse  Coupling  Coefficient  for  XeF  to  an  Aluminum 
Surface  as  a  Function  of  Incident  Fluence.  Energies 
of  0. 8-1.4  J  were  used. 

•  -  P  =  760  torr,  Tp  =  1.0  usee; 

o  -  P  =  75  torr,  Tp  =  1.0  ysec. 
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